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Abstract 

We carry out Monte Carlo simulations of the uniformly frustrated three 
dimensional XY model, as a model for vortex line fluctuations in high tem- 
perature superconductors in an applied magnetic field. We show, comparing 
systems of different size, that there are two distinct phase transitions. At a 
lower Tc_L, the vortex lattice melts, and coherence is lost in planes perpen- 
dicular to the magnetic field. At a higher Tcz, a vortex tangle percolates 
throughout the system, and coherence is lost parallel to the magnetic field. 
Cooling below Tcz, high energy barriers for vortex line cutting lead to an 

entangled glassy state. 
74.60.Ge, 64.60-i, 74.40.+k 
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I. INTRODUCTION 

In several earlier wor we have introduced the three dimensional uniformly frustrated 
XY model!, as a phenomenological model for studying phase transitions, and the effects 
of vortex line fluctuations, in the mixed state of high temperature superconductors in a 
uniform applied magnetic field. This model applies in the strongly type-II limit where 
the magnetic penetration length is much greater than the average vortex line separation, 
X ^ a^, and the magnetic induction inside the superconductor is approximately uniform. 
In our most recent work!, henceforth referred to as (I), we found evidence that the system 
undergoes two distinct phase transitions upon heating. First, at the lower Tc±, the vortex 
line lattice melts, destroying superconducting phase coherence in directions perpendicular 
to the applied magnetic field; coherence parallel to the field however remains. Then, at the 
higher T^^, coherence parallel to the field is lost as well. 

In the present work we extend the results of (I) in several major directions: (i) By 
presenting detailed studies of the system behavior as the system size is varied in the directions 
parallel and perpendicular to the magnetic field, we show clearly that the two distinct 
transitions found in (I) are not artifacts of finite size effects, (ii) We show that the upper 
transition T^z can be viewed as a percolation-like transition, where the vortex lines become so 
completely interconnected through mutual intersections, that one may trace out a connected 
path of vortex line segments which travels completely around the system in the direction 
perpendicular to the applied magnetic field, {in) We find, in contrast to our earlier results, 
that below Tcz the energy barrier for vortex line cutting grows so large, that cuttings are 
frozen out on the time scale of our simulation, and the system can cool into an entangled 
glassy state as in the "polymer glass" picture of Obukhov and Rubinsteini. 

The remainder of our paper is organized as follows. In Section II we outline our model 
and Monte Carlo method. In Section III we present the results of our simulations. Section 
IIIA gives results for the helicity modulus, which measures superconducting phase coherence. 
Section IIIB gives results for the average length of vortex lines due to thermal fluctuations. 
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Section IIIC analyses the entanglement of the vortex lines by considering the winding of the 
field induced vortex lines about the direction of the magnetic field. Finite size dependencies 
are investigated. Section HID gives results concerning the distribution of thermally excited 
closed vortex rings. Section HIE discusses the "2d boson" analogue to vortex line fluctu- 
ations, as applied to our model. In Section IV we summarize our results and discuss the 
possible connection to recent experiments. 

II. MODEL 

The model that we study is given by the Hamiltoniani 

n[e,] = JoY.v{e,-e,-A^) (i) 

where 9i is the phase of the superconducting wavefunction at site i of a three dimensional 
cubic numerical mesh, the sum is over all nearest neighbor bonds of this mesh, 

A„ = H^fA.dl (2) 

is proportional to the integral of the fixed magnetic vector potential A across bond {ij) 
($0 = hc/2e is the flux quantum), 

1 



\/(a) = -(T/Jo)ln exp 



Jo(a - 27rm)VT 



(3) 



is the Villain interaction between neighboring sites, and the coupling is 



where we identify the vortex core radius with the lattice constant a of our numerical mesh. 
For our numerical studies we take an isotropic uniform constant Jq, although this could 
be varied if desired to model the effects of anisotropy or randomness. Periodic boundary 
conditions are chosen in all directions. Numerical meshes of various sizes x Lz are studied 
(the subscript "_L" will refer to the x and y directions, transverse to the applied magnetic 
field). 



The approximations which lead from the famihar Landau-Ginzburg free energy functional 
to the Hamiltonian (|I|), and their justifications in the X ^ limit, have been described 
in detail in (I). In the following, we study the specific case where the uniform magnetic 
induction B = VxAis/ = l/25 flux quantum per plaquette of the mesh, oriented in the 
z direction. The ground state is a square periodic lattice of straight magnetic field induced 
vortex lines with spacing a^/a = l/vT = 5, as shown in Fig. |l]. Henceforth, we cite all 
lengths in units of a ~ ^o, and energies in units of Jq. The total number of field induced 
lines is = fL\. 

Our Monte Carlo simulations are carried out using the standard Metropolis algorithm. 
Performing the simulation in terms of the phase variables ^j, we locate the vortex lines in 
any particular configuration by computing the net phase change around every plaquette of 
the mesh. We define as an intersection, or cutting, between two vortex lines whenever we 
find a unit cell of the mesh which has more than one vortex line entering and leaving. In 
such a situation, we randomly assign which exiting segment is connected to which entering 
segment, for the purpose of identifying the paths of these particular lines. Each of our data 
points is typically the result of 2, 000 sweeps to equilibrate, followed by 15, 000 sweeps to 
compute averages, where each sweep refers to one updating pass through the entire numerical 
mesh. Our calculations were carried out on a Sparc 10 workstation; for our largest system, 
25^ X 200, each temperature point took approximately two CPU days of computation. 

III. RESULTS 
A. Helicity Modulus 

To investigate phase coherence, we compute the helicity moduli, which give the stiffness 
of the system to twisted phase boundary conditions (see (I) for derivation). In Fig. |^ we 
show our results for the the helicity modulus T ^{T) perpendicular to the magnetic field, and 
T z{T) parallel to the magnetic field, for lattices of fixed size L±_ = 25, but varying Lz = 50, 
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100, and 200. In Fig. ^ we show Tj_(T) and Tz(T) for approximately equal = 24, 25, but 
varying Lj_ = 25 and 50. We see clearly two transitions, with T_|_ vanishing at T^x ~ 1.35, 
and Tz vanishing at Tcz — 2.6. Comparing the results from different L± and L^, finite 
size effects are generally seen to be small, hence we have clear evidence for three distinct 
thermodynamic states. The middle state is one in which superconducting phase coherence 
is destroyed in planes perpendicular to the magnetic field, but coherence is preserved in the 
direction parallel to the magnetic field. Comparing the results for heating versus cooling, 
we see only a small hysteresis in T^, however hysteresis in T_l increases with increasing Lz- 
We will see that this hysteresis in T_l is related to the entanglement of the vortex lines as 
they cool into a glassy state. 



B. Vortex Line Lengths 

As a first measure of the amount and nature of vortex line fluctuations, we consider the 
average density of vortex line segments in different directions. In the ground state, the total 
length of vortex lines is N^Lz along z, while zero along x and y. In Fig. ^ we plot the 
vortex line length density Al^ versus temperature, where Al^ is defined as the total length 
of vortex line segments due to fluctuations (ie. in excess over the ground state value) in 
direction ft, normalized by N^L^. A/^ = 1 represents an excess vortex line length equal to 
that of the straight field induced lines at T = 0. We see that A/^ ^ ^ Alz in the vortex line 
lattice phase below Tc_l, as well as for much of the vortex line liquid phase between Tci_ and 
Tcz- This indicates that in these regions, the dominant fluctuations are directed transverse 
fluctuations of the magnetic field induced vortex lines, as shown schematically in Fig. 
Near and above Tcz however, we find that Al^^y — Alz. We will see that this is due to the 
proliferation of closed vortex ring excitations as shown in Fig. |^c. Above Tcz, Al^ > 1, and 
the total vorticity is dominated by the contribution from fluctuations. We will see that this 
region is an interconnected tangle of vortex line segments, with no unambiguous separation 
between field induced lines and thermally excited rings. 
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C. Entanglement 



To consider the entanglement of the field induced vortex lines, we make use of the periodic 
boundary condition which is imposed along the direction of the magnetic field z. If {r±i{z)} 
are the positions in the xy plane where the field induced vortex lines intersect the plane at 
constant z, then the set of points {r_Lj(0)} must be identical to the set of points {r±i{L;^)}. 
If we view this periodic boundary condition along z as representing the circumference of a 
three dimensional torus, then the magnetic field induced vortex lines will divide into distinct 
connected groups, each of which makes a certain number of windings around the system in 
the z direction before closing back on itself. A group making a winding m would consist 
of the m lines ii,i2,...,im satisfying the condition, rj_j^(0) = r±i^{Lz), r±i^{0) = r±i.^{Lz), 

r^j^(O) = r±i^{L;,). For example, in Fig. H, we show a configuration with two lines of 
winding m = 1, two groups of lines with winding m = 2, and one group of lines with winding 
m = 3. Qualitatively, when a configuration contains only windings with m = 1, we say that 
it is "unentangled" . When a configuration contains many windings with large values of m, 
it is highly entangled. To characterize the degree of entanglement, we compute the average 
distribution n(m) of the total number of lines n which participate in windings of value m. 
X]m^(^) = ^v, the total number of field induced vortex lines. 

In Fig. 1^ we plot versus T for several different system sizes, the ratio R = n{l)/N^ of lines 
which make a winding of m = 1. R = 1 indicates a completely unentangled set of lines. We 
see that for all sizes upon heating from the ground state, ~ 1 stays constant until about 
T ~ 2.0 > Tc_L, then decreases to its high T limit at T^z- Upon cooling however, R starts to 
rise below T^z and saturates around Tci_ to a value i? < 1, dependent on system size. Only 
for our shortest system, Lz = 50, do we find disentanglement, ie. R = 1 upon cooling. For 
all larger Lz, the lines remain trapped in a non-equilibrium entangled state upon cooling (in 
(I) our system size was Lz = 24, hence we failed to see the entanglement below Tcz that we 
now find). The degree of this entanglement increases (ie. R decreases) with increasing L^. 
As this cooled state is not in equilibrium, it is unclear if the low temperature value of R may 
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vary with independent coolings, or if it may strongly depend on the rate of cooling. To test 
this, we have carried out four independent coolings of the system size = 25, = 100, 
and find the T ^ values of R = 0.64, 0.68, 0.75, and 0.79. The second of these coolings 
was carried out using twice the number of Monte Carlo sweeps per temperature as for the 
rest of our data. We similarly have carried out four independent coolings of the system size 
L± = 15, L, = 100, finding values of R{T ^ 0) = 0.44, 1.0, 0.41 and 0.38. The second run 
indicates that large fluctuations in R are possible. 

The strong hysteresis we find in R, which measures the global topology of the lines, 
should be contrasted with the absence of hysteresis in the line length densities AZ^ (see 
Fig. ^), which are a local measure of line fluctuations (the slight hysteresis in Al^^y which 
appears below T ~ 1.0 < Tc± is due to the extra geometrical line length needed to make a 
quenched entangled state, compared to a lattice of straight lines). This suggests that the 
hysteresis in R is due to the energy barrier for the cutting of vortex lines. As T decreases 
below Tcz, thermal activation over this energy barrier, which is necessary to disentangle 
the lines, becomes frozen out on the times scales of our simulation. As a measure of this 
energy barrier, we compute the average number of line cuttings Nc (unit cells with more 
than one line entering and leaving) present in the system. We then define the cutting length 
= N^Lz/Nc as the average distance in the z direction between two successive cuttings of 
a single line. We plot versus T in Fig. H. Above T^z we find — 1 ~ 2 indicating a heavily 
interconnected tangle of lines with much cutting. As T decreases below T^^, increases 
rapidly, becoming of the order of L^. The absence of any size dependence in comparing 
the system with = 50 (which disentangles upon cooling) versus Lz = 200 (which remains 
entangled upon cooling) suggests that much of the cutting which determines in the region 
below T ~ 2.0 may be due to the intersection between field induced lines and thermally 
excited closed vortex rings, rather than between two field induced lines; cuttings between 
field induced lines may only be occurring on even larger length scales. This picture we 
find of cooling into a non-equilibrium entangled state is therefore similar to the "polymer 
glass transition" originally proposed by Obukhov and RubinsteinS. Comparing the data 
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for Lz = 50 with Lz = 200 in Fig. 0, it is interesting to note that with respect to phase 
coherence, entanglement has a noticeable effect only on the helicity modulus T^; seems 
entirely unaffected. 

As a further measure of the process of vortex line entanglement, we now consider the 
complete distribution of line windings n(m). In Fig. ^ we show n(m), for a fixed system 
size of L± = 25, Lz = 200, for various temperatures. The results shown were obtained upon 
cooling the system. We find that for all T > Tcz — 2.6, the distribution is n{m) = 1 for all m, 
ie. a field induced vortex line selected at random is equally likely to belong to a winding of 
any value m. This result is consistent with the assumption that each vortex line i is equally 
likely to reconnect onto any other vortex line j, upon traversing the system in the z direction 
once, ie. r_|_j(0) = rj_j{Lz) is equally likely for any i and j. The most likely explanation 
for such behavior is that above Tcz the lines become so completely interconnected due to 
cuttings, that the global path of a given line is primarily determined by our algorithm which 
makes a random choice for the continuation of the line at each individual cutting. When 
each line has sufficient cuttings with its neighbors, the resulting line path our algorithm 
traces out is equally likely to meander anywhere throughout the system. This conclusion is 
supported by Fig. ^ where we see — 1 "~ 2 for T > Tcz- We will see further evidence for 
this later when we consider the distribution of thermally excited vortex rings. 

As T decreases below T^z in Fig. ^ we find a steady increase in n{m) at smaller m, 
compensated by a decrease in n(m) at the largest m. This is as one would expect when the 
connectivity of any pair of field induced vortex lines i and j, ie. r_Lj(0) = r_i_j{Lz), becomes 
dominated by the thermal transverse wandering of the lines as they pass through the system 
along z, rather than by line cuttings. As T decreases, the transverse wandering decreases, 
and the probability for near neighbor reconnections increases with the resulting increase in 
n{m) for small m. It is interesting to note however, that even for T moderately below T^z, 
there remains a wide region of intermediate m, where we continue to find n{m) ~ 1. For 
T = 2.6 — 2.3 we believe that, even though we are below Tcz, the energy barrier for line 
cutting is still sufficiently small compared to T, that our data represents true equilibrium 
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behavior (see < 10 in Fig. and the absence of hysteresis in R in Fig. for these values 
of T). As T is cooled below ~ 2.3, the system gets trapped in some random metastable 
non-equilibrium tangle. 



We now consider the finite size dependence of entanglement. In Fig. ^ we plot n{m) 
for systems of fixed but varying = 50, 100, and 200. Our data is for the fixed 
temperature T = 2.4, below Tcz yet still high enough that we are sampling equilibrium. 
We see that as increases, n{m) approaches the T > Tcz limit of unity. This may be 
understood as a result of the increased transverse wandering of lines as increases, thus 
decreasing the probability of neighboring pair reconnections. Considering the value of ^(1), 
we see that it decreases by a factor ~ 2 as increases from 100 to 200, consistent with 
a random walk like behavior for the vortex line transverse fluctuations. This leads one to 
expect that in the limit ^ oo, for fixed L±, the system will remain completely entangled 
at all T < Tcz- 



In Fig. [TT] we plot n(m) at T = 2.4 for systems of fixed = 100, but varying Lj_ = 15, 
20 and 25. These are systems with a total of A^^, = 9, 16, and 25 field induced vortex 
lines respectively. As the maximum winding is always rrimax = N^, the fall off of n{m) 
at large m occurs at different m ^ for the different L±. If we normalize the different 
curves in Fig. [ll] by (recall, Z^m'^l"^) = ^v), we find that R = n{l)/Ny is approaching 
a constant value as oo (see also Fig. |^). Thus the fraction of disentangled lines is 

approaching a well defined value. However we show the curves without this normalization to 
point up the wide intermediate region where we continue to find n(m) ~ 1 as increases. 
These observations suggest that for T < Tcz for fixed L^, as increases, the equilibrium 
probability distribution of windings n{m)/Ny approaches a limiting form, but the average 
value of m diverges. 
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D. Vortex Ring Excitations 



We now consider the proliferation of thermally excited closed vortex ring excitations as 
illustrated in Fig. |^. Defining q{p) as the total number of vortex rings with perimeter p, 



we plot in Fig. the log of q{p) versus 1/T, for p = 2, 40. Our data is for the system 
size Lz = 100, Lj_ = 25. We show the results obtained from cooling; comparison with data 
from heating shows no significant hysteresis. For T < Tcz the data falls along straight lines 
over several orders of magnitude, clearly indicating a thermally activated form. These lines 
intersect at roughly the same temperature, I/Tq ~ 0.3, thus suggesting the low temperature 
form 

g(p)~goe-^(P)(^/^-^/"^°). (5) 



In Fig. we plot the the value of E{p), extracted from the data of Fig. |T2|, versus p and 
find the linear dependence, 

E{p) = -l.U + ep, £ = 3.32. (6) 

Thus for T < Tcz, the number of rings q{p) is determined by the excitation energy to create 
the ring, and this energy scales linearly with the ring perimeter. As T decreases, large rings 
get exponentially suppressed. 

For T > Tcz, we see from Fig. |l^, that q{p) saturates to a constant value, and that rings 
on all length scales p are now present. As discussed above in connection with the winding 
distribution n{m), we believe that this saturation of q{p) is the result of a transition in which 
the vortex lines become so heavily interconnected through cuttings, that a connected vortex 
tangle percolates through the entire system. In this heavily interconnected limit, there is 
in general no unambiguous way to classify a given vortex line segment as belonging to a 
particular ring of size p, or even as belonging to a ring versus a field induced line. The 
distribution q{p) would then be dominated by the statistics of our line tracing algorithm 
which makes random choices at each line cutting, rather than by any energetics. 
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In Fig. we replot our data as q{p) versus p for several different T. We show only data 
for sizes p in which the finite size effects, comparing different L±_ = 15, 20, 25, are small (to 
determine these finite size effects, we compared the normalized ring densities q{p)/LzL']_). 
For low T < Tcz — 2.6, we see an exponential decay q{p) ~ exp{—e'p/T), consistent with the 
discussion above. Comparison with Eqs.(|^) and (|^) gives for the effective ring line tension, 
e', at low temperatures, 

e' = il- T/To)e. (7) 

For larger T > Tcz, we see a slower than exponential decay, which is well fit by an algebraic 
power law, q{p) ~ p~^, x ~ 2.56, 2.69, 3.07, 3.75 for T = 5.0, 3.5, 3.0, 2.8 respectively. 
The cross over from exponential to algebraic decay occurs near T^z- The transition at Tcz 
can therefore be described as the vanishing of the ring line tension e' as T increases to Tcz- 
This picture has some similarities with proposed vortex ring unbinding theories of the phase 
transition in the ordinary three dimensional XY model@. 

E. 2d Boson Analogy 

As a final indication that Tcz is a vortex percolation-like transition, we compute a quantity 
motivated by Nelson's analogyil between the field induced vortex lines of a superconductor, 
and the imaginary time world lines of two dimensional bosons. According to this analogy, 
the 2d boson superfluid density ps is non-zero only when superconducting coherence parallel 
to the applied magnetic field is lost&i. A convenient expression for ps has been given 
by Cerpeley and Pollack0 in terms of the "winding number" W of boson world lines, 
Ps = mTbosoniW^) /2h^ where Tboson is the temperature of the boson system. The mappingll 
to the superconductor problem is given by: h /Tboson ^ Lz, Ti ^ Tguper, m ^ ei ~ ttJo the 
single vortex line tension. Hence ps ~ {W'^)/Lz. The winding number is defined^ in terms 
of the boson world lines, or equivalently in terms of the magnetic field induced vortex lines 
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1 

W = — E[r±.(^.)-rx.(0)]. (8) 

i=i 

W measures the net "winding" of the hnes about the system in the xy plane (W should 
not be confused with our earlier distribution n{m) which measures winding of lines about 
the z direction). Since the periodic boundary condition along z implies that the set of 
points {r_|_j(0)} is equivalent to the set of points {r^i{Lz)}, W can be non-zero only if 
periodic boundary conditions also exist in the x and y directions. If we assume that the 
only vortex lines present in the system are the magnetic field induced lines, then W is 
just equal to the net vorticity in the directions perpendicular to the magnetic field, or 
equivalently the perpendicular part of the q = Fourier transform of the vortex density 
n{r^,z) = (l/27r)V x V^, 

W = -^n^=o = J- [ng=o - z(z ■ ng=o)] • (9) 

Note a crucial difference between 11^=0 and the line densities Al^ we defined earlier: Al^ 
measures the total length of vortex line segments, independent of the direction of the vor- 
ticity; ng=o nieasures net vorticity, ie. two line segments oriented in opposite directions will 
cancel in their contribution to 11^=0 • 

The Hamiltonian of our system Eq. (|I|) can be expressedli in terms of the vortex density 

as, 

27^2 J 

'^[^q] = l—FT " /z<^9,o) ■ (ng - fz6g^o)Gg (10) 

J^z-Li± q 

where the interaction Gg ~ 1/g^ as g — >■ 0. To keep the total energy finite, we are thus 
rigorously constrained in our model to configurations where n^=o ~ 0- Hence as long as 
we assume that the only vortex lines present in the system are the field induced lines, we 
must have W = 0. If we now include the possibility of closed vortex ring excitations, the 
identification of Eq. (P) continues to be correct provided the rings remain of finite length p; 
the net vorticity of a finite ring always vanishes as the vorticity must always reverse direction 
in order for the ring to close back on itself. Only if we have rings so large (ie. infinite as 
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L±_ — i> oo) that they wind completely around the system in the x or y direction, making use 
of the periodic boundary conditions to close back on themselves without ever reversing the 
direction of their vorticity, will the identification between W and n^=o ® break down. 

A non-zero W computed as in Eq.(|^) is now possible, provided its contribution to n^^Q is 
exactly canceled by an oppositely oriented contribution to n^^Q from the infinite transverse 
ring. 

In Fig. |15| we plot (W'^)/Lz versus T, for system sizes = 100, L± = 15, 20, 25, and 
Lz = 50, L± = 25. We see that W'^ is only non-zero above Tcz — 2.6. Thus only above Tcz 
do we find vortex rings that travel completely around the system in the direction transverse 
to the applied magnetic field. This is only possible once the vortex tangle, of interconnected 
magnetic field induced lines and thermally excited rings, percolates throughout the entire 
system. Note that our results for {W'^)/Lz show some difficulties with the interpretation of 
this quantity as a 2d boson superfiuid density. Comparing sizes Lj_ = 25, Lz = 50, 100, we 
see no change in {W'^)/Lz, even though different Lz correspond to different temperatures 
Tboson/^ = 1/Lz in the 2c? boson problem. For fixed Lz = 100, and increasing we see 
a steady decrease in {W'^)/Lz towards zero, in contrast to expectations that the 2d boson 
Ps should approach a finite constant. We do not fully understand these size dependencies. 
It has been suggestedi that the results of Ceperley and Pollack for ps may not apply in 
the limit of a long range gauge interaction between 2d bosons, such as is the case in our 
superconductor problem. Nevertheless our results continue to support the view that T^z is 
a vortex percolation transition. 



IV. DISCUSSION 

Although we have not tried in this work to model a particular high Tc copper-oxide 
superconductor, it is worth indicating in what cases our results may qualitatively describe 
behavior in these materials. Our approximation of a uniform magnetic induction inside 
the material (formally equivalent to A ^ oo) means our model should apply only in the 
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limit where both the apphed magnetic field is large enougha that A ^ ~ y^o/B, and 
the Josephs™, oouphug betweeu the CuO planes do,m„ates0 over the magnetic couphng, 

Ac < A^f,/(i (where Ac, Xab are the magnetic penetration lengths perpendicular to and within 
the CuO planes respectively, and d is the separation between the CuO planes). Behavior 
in the high Tc materials has been further characterized with reference to a critical magnetic 
fieldM, Bcr — $0'^ab/-^c'^^- ^or B < B^r, vortices within different CuO planes may be 
thought of as correlated strings, an anisotropic Landau-Ginzburg description is adequate, 
and the melting of the vortex line lattice is "three dimensional". For B > B^r, vortices 
within different CuO planes are weakly coupled, the layered Lawrence-Doniach mo del0 IS 
more appropriate, and melting is "quasi-two dimensional." Since in our simulation we have 
taken Jq ~ 1/A^ constant in all directions, and the spacing between vortices is a^/a = 5, 
if we identify the lattice constant of our numerical mesh a with the spacing between CuO 
planes d, we have ^o/al = B < B^r and our results apply in the region where melting is 
three dimensionallii. 

One of the primary results of this paper has been to substantiate the existence of two 
distinct phase transitions in our model. At the lower transition Tc±, the vortex line lattice 
melts upon heating (see (I)), — > 0, and superconducting coherence is lost in the planes 
perpendicular to the applied magnetic field. At the upper transition Tcz, a vortex tangle 
percolates completely through the system in the directions transverse to the applied magnetic 
field, ^ 0, and superconducting coherence is lost in the direction parallel to the magnetic 
field. Identifying the loss of superconducting coherence with the onset of linear electric 
resistivity, we therefore would expect the following experimental consequences: As T is 
decreased, the linear resistivity for currents applied parallel to the magnetic field will vanish 
below Tcz- However linear resistivity for currents applied perpendicular to the magnetic 
field will continue to remain finite below Tcz until a lower Tc± is reached. This result is 
in agreement with predictions by Feigel'man and co-worker£i, as well as by Glazman and 
Koshelev0, for B < Brr- 



14 



Recent experiment sl3 by Steel, White and Graybeal on synthetic MoGe/Ge multilayers 
appear to show precisely such behavior. In these experiments, in which the magnetic field is 
applied perpendicular to the layers, the authors observe a well defined temperature "To" at 
which the resistivity parallel to the magnetic field shows a dramatic drop, accompanied by 
the onset of substantial nonlinearities in the I — V characteristics. This suggests a transition 
where the linear resistivity in this direction vanishes. The resistivity perpendicular to the 
magnetic field shows a kink at To, however continues to remain linear for temperatures 
T < To- Such behavior is consistent with that of our middle phase Tci_ < T < Tcz, if we 
identify the experimental To with our Tcz- These experiments however appear to be in the 
region B > B^r, so the direct application of our results remains unclear. 

A second important result of our paper has been the observation that upon cooling below 
Tcz, lines can get trapped in a disordered entangled state where vortex line cutting is frozen 
out except on long time scales. Here the disorder is purely topological in nature^ and not 
due to any random impurities. The importance of such entanglement on transport prop- 
erties determined by vortex line diffusion, has been stressed by Nelson and co-worker£0, 
particularly with regard to pinning by large scale impurities. 

Finally, we have identified the upper transition of our model, Tcz, as the temperature at 
which an interconnected tangle of wandering vortex lines and thermally excited vortex rings 
percolates through the system. Our analysis of the vortex ring distribution suggestes that 
this is the temperature at which the effective vortex line tension vanishes. Since a vortex 
line may qualitatively be viewed as a one dimensional "interface" between different ground 
states of the 3d XY model, the two transitions of our model might be viewed in analogy to 
the behavior of interfaces in the 3d Ising model. Our lower melting transition Tc± might be 
viewed like a "roughening" transition. Below Tc± lines remain straight as ^ oo, and are 
periodically ordered. Above Tc± line wandering increases with increasing Lz, and lines are 
disordered in the plane; however lines retain a finite line tension and so remain well defined 
fiuctuating objects. Our upper transition Tcz niight be viewed as the "bulk" transition where 
the effective line tension vanishes, and detached "bubbles," ie. vortex rings, proliferate. A 
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similar picture is implied in work by Bulaevskii et a/. Ill, and by Glazman and Koshelev0. 

In our discussion of the 2d boson analogy in section HIE, we derived the important con- 
sequence that ps = for all T below the vortex percolation transition, from the observation 
that the vortex line interaction of our model was Gq ~ l/q"^, and hence energy conserva- 
tion strictly requires W ~ ii^=o = 0- This is a direct consequence of our approximation 
A oo. For finite A, the interactioni is Gq ~ + A^^), and now fluctuations with 

finite n^^Q > are energetically possible in the vortex line hquid phase. This could imply 
finite ps, and vanishing superconducting coherence along the direction of the magnetic field. 
Recent work by Tao and Teiteli however shows that, for finite A, superconducting coherence 
along the magnetic field should still persist in a hexatic vortex line liquid state0, which 
might exist intermediate to the vortex line lattice, and normal vortex line liquid states. In 
this case, it remains to be seen if the hexatic to normal line liquid transition coincides with 
the vortex percolation transition of our model, or if the vortex percolation remains a sharp 
thermodynamic transition at all. 
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FIGURES 

FIG. 1. Ground state vortex line lattice for a magnetic induction of / = B^q/^q = 1/25 flux 
quantum per unit cell of the numerical mesh. The view is along the direction of B and (+) locates 
the positions of the straight vortex lines. 

FIG. 2. Helicity modulus along the direction of B, and perpendicular to B, for lattice 
sizes L± = 25 and varying = 50, 100, and 200. Both heating and cooling are shown. The 
vanishing of '^z,± indicates two separate transitions. No significant finite size effects are seen. 

FIG. 3. Hehcity modulus along the direction of B, and perpendicular to B, for lattice 
sizes Lz = 25, L± = 25 and = 24, L± = 50. Both heating and cooling are shown. No significant 
finite size effects are seen. 

FIG. 4. Line length densities Al^, measuring absolute value of total vortex line lengths in 
direction /t, normalized by total length in ground state NyL^. For AZ^, the ground state line length 
has been subtracted, in order to show only excess length due to fluctuations. Both heating and 
cooling are shown for a fixed lattice size L± = 25, = 200. The solid horizontal line indicates the 
total normalized length in the ground state; the solid vertical lines mark the transition temperatures 
as obtained from the vanishing of 

FIG. 5. Schematic of possible vortex fluctuations, (a) shows a directed fluctuation of a field 
induced vortex line; the line pierces each plane of constant z only once. (6) shows a field induced 
vortex line with an overhang, (c) shows a closed vortex ring excitation. 

FIG. 6. Schematic of possible reconnections of field induced vortex lines, under application of 
the periodic boundary condition in the z direction. This example shows two lines of winding m = 1, 
two groups of lines with winding m = 2, and one group of lines with winding m = 3. Solid, dashed, 
and dotted lines are use to distinguish the different lines within a particular winding group. 
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FIG. 7. Fraction of field induced vortex lines which axe unentangled, R = n(l)/Ar„, for various 
system sizes. Note the strong hysteresis between cooling and heating. Entanglement increases as 
Lz increases. 

FIG. 8. Distance along z between two successive cuttings of a single field induced vortex 
line, '"^ 2 — 1 for r > Tcz indicates a heavily interconnected vortex tangle. 

FIG. 9. Distribution of windings n{m) that field induced vortex lines make in traveling around 
the system along the z direction. Several different temperatures are shown for the fixed system 
size = 15, Lz = 200. For T > T^^ ~ 2.6, we find n(m) = 1. 

FIG. 10. Distribution of windings n(m) for fixed T = 2.4 < Tcz, for system sizes L± = 25 and 
Lz = 50, 100, 200. As Lz increases, n(m) — 1, ie. entanglement increases. 

FIG. 11. Distribution of windings ra(m) for fixed T = 2.4 < Tcz, for system sizes Lz = 100 and 
L± = 15, 20, 25. Distribution remains flat, n(m) ~ 1, for wide region of intermediate m as 
increases. 

FIG. 12. Distribution of thermally excited closed vortex rings of perimeter p, versus l/T, for 
fixed system size L± = 25, Lz = 100. Straight solid lines for T < Tcz show thermally activated 
behavior. 

FIG. 13. Energy barrier E(p) for vortex rings of perimeter p, as extracted from the T < Tcz 
data of Fig. 12. E{p) scales linearly with p. 

FIG. 14. Distribution of thermally excited closed vortex rings of perimeter p, versus p. Several 
different temperatures are shown for the fixed system size L± = 25, Lz = 100. For T < Tcz — 2.6 
solid lines are the best fit to an exponential decay. For T > Tcz solid lines are the best fit to an 
algebraic decay. 

FIG. 15. Winding number {W'^)/Lz versus temperature for system sizes Lz = 100, L± = 15, 
20, 25, and Lz = 50, L± = 25. >() only for T > ~ 2.6. 
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